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ABSTRACT 

We explore the mult i- messenger signatures of encounters between two neutron stars 
(ns^) and between a neutron star and a stellar- mass black hole (nsbh). We focus on the 
differences between gravitational wave driven binary mergers and dynamical collisions 
that occur, for example, in globular clusters. Our discussion is based on Newtonian 
hydrodynamics simulations that incorporate a nuclear equation of state and a multi- 
flavour neutrino treatment. For both types of encounters we compare the gravitational 
wave and neutrino emission properties. We also calculate the rates at which nearly 
unbound mass is delivered back to the central remnant in a ballistic-fallback-plus- 
viscous-disk model and we analyze the properties of the dynamically ejected matter. 
Last but not least we address the electromagnetic transients that accompany each 
type of encounter. 

We find that dynamical collisions are at least as promising as binary mergers for pro- 
ducing (short) gamma-ray bursts, but they also share the same possible caveats in 
terms of baryonic pollution. All encounter remnants produce peak neutrino luminosi- 
ties of at least ^ 10^^ erg/s, some of the collision cases exceed this value by more than 
an order of magnitude. The canonical ns^ merger case ejects more than 1% of a solar 
mass of extremely neutron-rich (Y^ ~ 0.03) material, an amount that is consistent 
with double neutron star mergers being a major source of r-process in the galaxy, 
nsbh collisions eject very large amounts of matter (~ 0.15 M©) which seriously con- 
strains their admissible occurrence rates. The compact object collision rate (sum of 
ns^ and nsbh) must therefore be less, likely much less, than 10% of the ns^ merger 
rate. The radioactively decaying ejecta produce optical-UV "macronova" which, for 
the canonical merger case, peak after ~ 0.4 days with a luminosity of ~ 5 x 10^^ erg/s. 
ns^ (nsbh) collisions reach up to 2 (4) times larger peak luminosities. The dynamic 
ejecta deposit a kinetic energy comparable to a supernova in the ambient medium. 
The canonical merger case releases approximately 2 x 10^^ erg, the most extreme (but 
likely rare) cases deposit kinetic energies of up to 10^^ erg. The deceleration of this 
mildly relativist ic material by the ambient medium produces long lasting radio flares. 
A canonical ns^ merger at the detection horizon of advanced LIGO/Virgo produces a 
radio flare that peaks on a time scale of one year with a flux of ~0.1 mJy at 1.4 GHz. 
Collisions eject more material at higher velocities and therefore produce brighter and 
longer lasting flares. 

Key words: black hole physics - gravitational waves - neutrinos - nuclear reac- 
tions, nucleosynthesis, abundances - radiation mechanisms: non-thermal — gamma-ray 
bursts 
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The encounter of a neutron star (ns) with another neutron 
star or with a stellar mass black hole (bh) is a fascinating 
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Figure 1. Initial density and Ye profiles (in hydrostatic and /3-equilibrium) of the neutron stars used in this study (1.3 and 1.4 M©, 
nuclear equation of state of Shen et al. ( |Shen et al.|1998a| ). 



events that involve a large variety of different physical pro- 
cesses, see 



Faber I p009\ , [Duez] ([2010t , [Shibata k TanigucE 
( 20111 , |Rosswog| V2011') and'F aber Rasio| ( |2012| ) for recent 



reviews on various aspects of this topic and for a guide to 
the current literature. 

The detection of the gravitational waves (GWs) emitted dur- 
ing the inspiral of ns^ and nsbh binary mergers are prime 
targets of ground-based gravitational wave detectors such as 



LIGO, Virgo and GE06QQ (Willkeet al.|2QQ6; Acernese et al. 



2008| |Grote||2QQ8| |Smith||2009t . In their advanced states, 
the interferometers target to detect the signals of ns^ coa- 
lescences out to hundreds of megaparsecs, corresponding to 
redshifts of z ~ 0.1. Once detected, the gravitational waves 
offer the exciting possibility of filling gaps in our understand- 
ing of the neutron star equation of state in the high density, 
low temperature regime that is experimentally hardly acces- 
sible. 

Compact binaries have for a long time been the prime can- 
didates for the central engines of (short) gamma-ray bursts 
(GRBs; [Eichler et al.| (|1989|); |Paczyhski| (|1991|); [Narayan 



et al. (1992)) and this hypothesis has survived being con- 
fronted with a wealth of observational results in recent years. 
Several challenges remain, however, and the case is far from 



being closed (,Piran 


2005 


jNakar 1^2007 Lee & Ramirez-Ruiz 


2007 Gehrels et al. 


2009 


). 



They are also promising sources for the heaviest elements 
in the Universe that are formed via rapid neutron cap- 
ture ([Lattimer fc Sch ramm 19 741 1 1976 1 fEichler et al.||1989| 
Freiburghaus et al.|1999,:,Rosswog et al.|1999| [Roberts et al.| 



2011 1 Goriely et al.| 2011). The textbook r-process source, 
core-collapse supernovae, have been found to be seriously 
challenged in providing the physical conditions (high en- 
tropy, low electron fraction together with rapid expansion) 
that are required to produce the heavy (A > 90) r-process 
elements ( Roberts et al.pOlO Fischer et al.|[20To] [Arcones 



Janka|[20TT| H The main r-process nucleosynthesis con- 
tenders are compact binary mergers which release neutron- 
rich matter in at least three ways. Apart from the matter 
that is ejected dynamically via gravitational torques, there 
is an additional contribution due to neutrino-driven winds 
( Dessart et al.|2009 ) and from the late-time dissolution of ac- 



cretion disks (Chen Beloborodov|2007||Beloborodov|2008 



[Metzger e t al.||2008 ). While the initial starting point is the 
same, cold neutron star matter in /^-equilibrium, the three 
channels differ in the amounts of released matter, in their en- 
tropies, expansion time scales and electron fractions. There- 
fore they might possibly produce different nucleosynthetic 
signatures. 

The ejecta are responsible for two types of electromagnetic 
transients: the decompressed neutron star matter is sub- 
ject to radioactive decays (Li & Paczyhski 1998| Rosswog 
|2005| [Metzger et al.|2010| [Roberts et al.|2011t and at some 
stage the ejecta dissipate their kinetic energy in the am- 
bient medium. The former is expected to produce an op- 
tical display not too different from a supernova, but much 
shorter, sometimes referred to as "macronova", the latter 
has been shown to produce detectable, isotropic radio emis- 
sion that peaks near one gigahertz and persists on a de- 
tectable (sub-milliJansky) level for weeks out to a distance 
oi z^O.l dNakar k Piran|2011[ ). Especially if the true GW 
detection rates should be near the lower end of the predic- 
tions, additional electromagnetic signatures would be crucial 
to confirm marginal gravitational wave detections and they 
would therefore enhance the effective detector sensitivities 



dKochanek Piran| 19931 [Hughes ^ Holz|2003[ [Dalai et al. 
2006|[Arun et al.[[2009[ ) 



^ A possible exception may be magnetorotationally driven super- 
nova jets where interestingly low electron fraction values seem to 
be reachable ( [Winteler et al.[[2012| . It remains to be explored, 
however, how robust this scenario is with respect to the stellar 
parameters and with respect to its nucleosynthetic yields. 
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Recently, dynamical collisions between compact objects as 
they may occur in the core of globular clusters have been 
studied by several authors ( [Kocsis et al. 2006 O'Leary 



et al. 



2009 



Lee et al. 2010| Kocsis & Levin 



et al. (2010) concluded that collisions could produce GRBs 
at a detectable rate. To demonstrate the viability as GRB 
engines they performed hydrodynamic simulations of com- 
pact object collisions, though without use of detailed mi- 
crophysics such as a nuclear equation of state or neutrino 
emission. 

The main questions that we want to address in this paper 



a) In which ways do the remnants and signatures of dy- 
namical collisions differ from those of binary mergers? 

b) To which extent can their rate be constrained by nucle- 
osynthetic yields? 

c) How different are electromagnetic transients following 
mergers and collisions? 

To address these issues we performed a sizeable set of New- 
tonian hydrodynamics simulations which include a nuclear 
equation of state and an opacity-dependent neutrino cool- 
ing scheme. These simulations are subsequently explored to 
predict the electromagnetic signatures of compact binary 
encounters. Clearly, the investigated systems are relativistic 
and ultimately General Relativity should be applied. It is, 
however, not gravity alone that shapes the dynamics and 
observable signatures of these encounters. Instead, also the 
remaining fundamental interactions contribute their share: 
a) the strong interaction via the nuclear equation of state, b) 
the weak interaction since it determines the neutrino emis- 
sion rates and thus the evolution of the electron fraction Ye 
and c) the electromagnetic interaction which is, for exam- 
ple, responsible for producing radio-flares once ejected ma- 
terial dissipates its kinetic energy in the ambient medium. 
Given this complexity, we consider Newtonian gravity as a 
tolerable approximation for the time being. The presented 
simulations are meant to serve as benchmarks for future sim- 
ulations that may include general relativity and the relevant 
microphysics. 

We discuss the viability of dynamical collisions as GRB cen- 
tral engines, in particular the properties of remnant disks, 
their neutrino emission and the prospects for magnetic field 
amplification. We further calculate the mass and the return 
time scales of fallback and we present in detail the properties 
of dynamical eject a as a basis for subsequent nucleosynthesis 
calculations. Based on these findings we discuss the question 
"What is the electromagnetic signature of a ns^ and a nsbh 
encounter?" 

The paper is structured as follows. In Sec. [2] we briefly sum- 
marize the ingredients of our simulations both in terms of 
physics and numerical methods. In Sec. [3] we describe the 
main findings. Sec. |3.1| discusses the dynamics and its im- 
print on gravitational waves. Sec. |3.2| explores the neutrino 
signal and Sec. |3.3| the differences between merger and col- 
lision fallback. In Sec. |3.4| we discuss the properties of the 
dynamically ejected material such as mass, electron fraction 
and velocity structure. They all shape the electromagnetic 
display, which is addressed in Sec. |3.5| Our results are sum- 
marized and discussed in Sec. [U 



2 SIMULATIONS 

The simulations of this paper make use of the Smooth Par- 



ticle Hydrodynamics (SPH) method, see Monaghan (2005) 
and Rosswog] (j2009| for recent reviews. Our code is an up- 
dated version of the one that was used in earlier studies 



(Rosswog & Davies 2002; Rosswog & Liebendorfer 



Rosswog et al.||2003| |Rosswog||2005 ). We solve the follow- 



2003 



ing evolution equations for each particle a 



dVg 

dt 

dUa 

"df 

dYe,t 
dt 



^^b{^ + ^+Uab] VaWab + /a,g + faM 



pa Pb 



rUb ( ^ + ^Uab I Vab • VaVl/^a 



dUa 



dt 



Apc Yn — AeC 



(3) 



the mass density is calculated by summing up contributions 
from neighboring particles 



Pa = y^^rribWab- 

b 



(4) 



Here rub is the (constant) mass of particle h and Wab 
W(\fa — rb\,hab) denotes the cubic spline kernel (Mon- 



|aghan|p"985 ) evaluated with the average smoothing length 
hab = {ha + hb)/2. fa,g IS the additional acceleration due 
to self-gravity that we evaluate using a binary tree ( Benz| 
|et al.||l990| and fa ,GW results from gravitational wave 
back-reaction. The relative particle velocity is denoted as 
Vab = Va — Vb. To producc entropy in shocks artifical dissi- 
pation is included via the tensor Hab. It has the standard 
form Monaghan] ( |1992 ) but particular care has been taken 
to avoid possible artifacts due to artificial viscosity, this has 
been outlined in detail in Rosswog et al. (2008). The quan- 



tity Ua denotes the specific internal energy of particle a the 
evolution of which is determined by PdV-work and viscous 
heating (summation term) and the energy loss to neutrinos, 



dt 



The quantities 



Apc 



oeff 



^eff 



and Aec 



(5) 



are the electron and positron capture rates per neu- 
tron/proton. R^c/pc effective neutrino number 
emission rates and the quantities r/np/ripn reduce in the non- 



degenerate limit to the number densities Un and Up ( Bruenn 



1985), for a more detailed account on the neutrino treat- 



ment we refer to the original paper ( Rosswog & Liebendorfer 



2003). The pressure at a particle b, Pb{pb,Tb,Ye,b), is calcu- 



lated using the Shen et al. equation of state (EOS) (Shen 
et al. 1998a|b ) extended to lower densities as described in 



Rosswog k Davies (2002) 



Compact binary systems are driven towards coalescence via 
the emission of gravitational waves. The corresponding ra- 
diation reaction forces for a slow-motion, weak field source 
can be calculated as the gradient of a radiation reaction 
potential which contains the fifth time derivatives of the re- 
duced quadrupole moments ( Burke|1971 ). Simple backreac- 
tion prescriptions usually rely on reducing the order of the 
time derivatives by averaging over several orbital periods. 
This procedure is well justified during the secular inspiral 
stages of a compact binary merger, but it has no well-defined 
meaning in the case of a parabolic encounter. For this reason 
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we ignore the exerted backreaction for the collision cases of 
this first study, fo r the merger cases w e use a simple point 



mass description (Davies et al. 199?) for /a,GW- The pos- 



sible impact of this technical shortcoming will be discussed 
below. 

As a diagnostics of the dynamical evolution we will use be- 
low the gravitational wave amplitudes as seen by an ob- 
server located along the binary rotation axis. Consistent 
with the Newtonian treatment of gravity we extract gravi- 
tational waves in quadrupole approximation 



+ 

TT 



IG 
2G 



(r 



(6) 
(7) 



where P-^ is the reduced quadrupole moment tensor evalu- 
ated at retarded times, 



(8) 



The needed time derivatives can be obtained by straight 
forward differentiation of Eq.(|8| so that the amplitudes can 
be calculated as sums involving particle masses, positions, 
velocities and forces. 

In the presented simulations we restrict our collision study 
to parabolic orbits. The strength of such an encounter is 
parametrized by the parameter 



/3 



rp 



where the Ri are the neutron star/Schwarzschild radii of 
the involved objects and rp is the pericenter distance. Thus 
/3 = 1 corresponds to a grazing impact, stronger (weaker) 
impacts have larger (smaller) values. Since the collision rates 
are proportional to the pericenter distance rp, we consider 
run A with /3 = 1 as the most likely ns^ collision case. Col- 
lisions with pericenter distances /3 < 1 can still form tidal 
capture binaries ( Fabian et al.||1975| |Lee et al.||2010| ) which 
can lead to final collisions after a sequence of pericenter 
passages. Collisions with /3 < 1, however, become compu- 
tationally increasingly cumbersome due to the discrepancy 
between orbital and internal dynamical time scales. There- 
fore we only consider collisions with /3 ^ 1. To keep the 
explored parameter space under controle, all investigated 
nsbh collisions possess a fixed impact strength of /3 = 1 and 
we only vary the black hole mass. 

It has long been known that the ns mass distribution pos- 
sesses a narrow peak near 1.35 Mq (|Thorsett & Chakrabarti 
19991 . Recent 
around 1.5-1.7 



studies 

Mq 



find an additional broader peak 



(Kiziltan et al. 2010 Valentim et al 



2011 ) for neutron stars with white dwarf companions. There 
may be an additional low-mass peak near 1.25 M0 produced 
by electron capture supernovae (Podsiadlowski et al. 2004; 



van den Heuvel||2004| [Schwab et al.||2010| . The niass dis- 



tributions for neutron stars of d ifferent evolutiona ry paths 
have recently been discussed in 



Ozel et al. 



(2012). In our 



simulations we restrict ourselves conservatively to masses 
near the 1.33 Mq peak that the latter authors find for dou- 
ble neutron star systems. 

In ns^ collision cases we use masses of 1.3 and 1.4 and 
explore the dependence on the impact strength parameter /S. 
For nsbh collisions we use a neutron star of 1.3 Mq and we 



investigate the dependence on the black hole mass (ruhh = 3, 
5 and 10 M0) while keeping the impact parameter at /3 = 1. 
We consider run H with mi = 1.3 M©, 777,2 = 1.4 M©, 
i.e. q ^ " 



0.923 an d negligible spins ( [Bildsten Cutler [1992 

Kochanek||1992 ) as the generic merger case and we will use 



it frequently as a reference point to compare the other cases 
against. As a somewhat academic case, we explore an ini- 
tially tidally locked binary neutron star system in run G. In 
runs I and J we briefiy touch upon neutron star black hole 
mergers where the black hole is treated as a Newtonian point 
mass with an absorbing boundary at the Schwarzschild ra- 
dius. 

The initial neutron stars are constructed from spheri- 
cally symmetric, zero-temperature, /^-equilibrium profiles, 
see Fig. [l] The SPH particles are placed in a close packed, 
hexagonal lattice configuration and they are subsequently 
relaxed so that they can find their true numerical equilib- 
rium state, see Sect. 3.1 of [Rosswog Price (2007). Some 



of our simulations have been run up to 0.5 s, more than an 
order of magnitude longer than existing simulations on this 
topic. 

All performed simulations are summarized in Tab. [l] 



3 RESULTS 

3.1 Encounter dynamics and gravitational wave 
emission 

To set the stage for later comparisons, we start with a brief 
description of the "standard" binary merger case, run H. In 
Fig. [2] we show a 3D rendering of its temperature distribu- 
tion. We only display matter below the orbital plane so that 
the temperatures and fiow structures inside the central rem- 
nant can be easily grasped. About one orbital period after 
contact two asymmetric spiral arms have formed (panel one 
and two), which evolve during the next ^15 milliseconds 
into a nearly axisymmetric torus (panel 4). When the stars 
come into contact a shear interface forms between them. 
Such Kelvin-Helmholtz unstable interfaces have long been 
known to emerge in neutron star mergers, see, for example, 
'Ruffert et ar|(|1996t; [Rosswog et ar|(|1999t;|Rasio &: Shapiro 



( 1999|. The resulting vortices have also been found to locally 
amplify pre-existing magnetic fields (Price Rosswog|2006 



,Anderson et al.,,2008 i Obergaulinger et al.||2010 ) and inside 



of them the (SPH particle) temperatures can temporarily 
reach values in excess of 60 MeV. The somewhat academic 
case of an initially tidally locked binary shows more pro- 
nounced tidal tails (due to larger angular momentum), but 
similar temperatures. Both double neutron star merger cases 
produce reasonably well-defined massive tori of 0.25 Mq in 
the irrotational "standard" case and 0.30 Mq for tidal lock- 
ing, see Tab. [2] 

The collision cases in contrast can suffer several close en- 
counters before finally merging into a single object and dur- 
ing these passages the neutron stars are efficiently tidally 
spun up. In the /3 = 1 case a single object only forms after 
the third close encounter, see Fig. [3] In the first, grazmg 
impact the stars' obital energy is used to spin up the stars 
to close to their breakup period, e.g. panel 2, and now they 
form an eccentric tidal capture binary. The next encounter 
near t = 8 ms is more central and again produces strong 
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Table 1. Overview over the performed simulations. The impact strength parameter /S is defined in Eq. (|9|. 



Run mi [Mo] m2 [M©] /3 A^sph [10^] tend [ms] objects/comment 



Colhsions 














A 


1.3 


1.4 


1 


2.7 


21.2 


ns-ns 


B 


1.3 


1.4 


2 


8.0 


9.0 


ns-ns 


C 


1.3 


1.4 


5 


2.7 


13.2 


ns-ns 


D 


1.3 


3.0 


1 


1.3 


127.5 


ns-bh 


E 


1.3 


5.0 


1 


1.3 


143.6 


ns-bh 


F 


1.3 


10.0 


1 


1.3 


540.3 


ns-bh 


Mergers 














G 


1.3 


1.4 


n.a. 


2.7 


20.3 


ns-ns, corot. 


H 


1.3 


1.4 


n.a. 


2.7 


19.1 


ns-ns, no spins 


I 


1.4 


5.0 


n.a. 


0.2 


138.7 


ns-bh, no spins 


J 


1.4 


10.0 


n.a. 


0.2 


139.3 


ns-bh, no spins 



Kelvin-Helmholtz vortices at the interface in which (SPH 
particle) temperatures locally exceed 80 MeV. The stars 
separate once more, with the 1.3 M0 star now transferring 
mass in a direct impact phase into the primary, see panel 4. 
The final encounter occurs around t ~ 12 ms, again forming 
a string of Kelvin-Helmholtz vortices (panel 5) and finally 
shedding mass from the secondary neutron star (panel 6). 
During the encounter the density never exceeds the initial 
value 3.6 x 10^^ gcm"^). 

The more central encounters form a single object after two 
{/3 = 2) and just one encounter {/3 = 5). In both cases 
strong shocks form in which the (SPH particle) temperatures 
reach values in excess of 80 MeV. In such shocks the neutron 
stars are substantially compressed, to values of ^ 4.52 x lO^'^ 
gcm~^(/3 = 2) and ^ 5.55 x 10^^ gcm~^(/3 = 5). In a super- 
position of rapid rotation and violent, stellar-radius ampli- 
tude oscillations, the central objects produce a multitude of 
interacting shocks in a string of mass shedding episodes, see 
Fig. [4] The oscillations are also imprinted on the neutrino 
signal, see below. 

For the neutron star black hole cases we also begin with 
binary mergers as a reference point. Since they have been 
explored in detail before ( Rosswog et al. 2004||Rosswog 2005 



Rosswog|2007bt , we restrict ourselves to a brief summary. In 
the case of a 1.4 Mq ns and a 5 Mq bh (run I) the neutron 
star starts transferring mass into the hole after 1.5 orbital 
periods. Consistent with our earlier studies this does not 
lead to the disruption of the neutron star on a dynamical 
time scale. Instead, self- gravity overcomes tidal forces again 
and the neutron star enters a long-lived phase of episodic 
mass transfer during which it transfers mass periodically 
towards the hole while shedding mass through its outer La- 
grange poinl|^ This phase continues for as many as 25 orbital 
revolutions before the neutron star is finally completely dis- 
rupted. The remnant at the end of the simulation (t = 138.7 
ms) consists of a "disk inside a disk" with a mass of 0.16 Mq 



^ Phases of stable mass transfer are not restricted to the case of 
Newtonian gravity. A stiff equation of state ( [Rosswog et al.|2004| ), 
small mass ratios and large bh spin parameters make systems par- 
ticularly prone to stable mass transfer, see ,Shibata Taniguchi, 
(20111 for a further discussion. 



for the inner, high density disk (p > 10^^ gcm""^, r < 120 
km) and 0.22 M© if also the outer disk {p > 10^ gcm~^, 
r < 700 km) is counted. The dynamics of the 1.4 M0(ns) - 
10 M0(bh) system proceeds in a similar manner, here after 
15 orbital revolutions the neutron star is finally disrupted 
and leaves a 0.20 disk together with a rapidly expanding 
one armed spiral structure. All the numerically determined 
mass trasnfer durations must be considered as robust lower 
limits on the true values (D an et al.|2011 ). 
For the neutron star black hole collision cases we only ex- 
plore the dependence on the black hole mass and keep the 
impact strength (/3 = 1) and neutron star mass (mns = 1.3 
M0) constant. During the first pericenter passage of the 
mbh = 3 Mq case, run D, the neutron star survives as 
a tidally spun up (close to break up, P ~ 0.95 ms) self- 
gravitating object, but sheds some of its mass in a tidal tail. 
When the neutron star passes the black hole after about 5 
ms for a second time another tidal tail is produced. Once 
more, the core of the neutron star survives as a gravitation- 
ally bound object. It is only completely disrupted during the 
third and final pericenter passage at t ~ 11 ms. After 127 ms 
the remnant consists of the bh with 3.98 M©, surrounded by 
a massive disk 0.15 M0, see Tab.|2]) which is externally 
fed by three spiral arms. Qualitatively, the mbh = 5 Mq case 
evolves in a similar manner, see Fig.[5l but now the neutron 
star core survives even the third passage. At the time when 
we have to stop the simulation (t = 144 ms), the core is, ac- 
cording to its radial velocity, unbound from the black hole. 
The neutron star core, however, is embedded into the debris 
gas and might therefore be further braked during its subse- 
quent evolution so that it will possibly fall back towards the 
bh. In the mbh = 10 M© case the neutron star is already 
completely disrupted during the second pericenter passage. 
Consistent with the findings of Lee et ah] ( |2010 ), all nsbh 
encounters have in common that they all leave behind a bh 
with a massive remnant disk (see Tab. [2| and one tidal tail 
per close encounter. 

The orbital dynamics is imprinted on the gravitational wave 
(GW) signal, for its calculation see Sec. [2] For the neutron 
star mergers (run G and H) the gravitational wave ampli- 
tudes /i+ (times the distance to the source d as measured 
in code units of 1.5 km) are shown in Fig. |6] upper left 
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Figure 2. 3D rendering of the temperature distribution in the standard neutron star merger case (1.3 and 1.4 Mq, no spin; run H). The 
upper half of the matter distribution has been "chopped off' to allow for a view into the stars. To enhance the contrast, the upper limit 
of the colourbar has been fixed to 20 MeV. In the various vortices that emerge due to Kelvin-Helmholtz instabilities peak temperatures 
in excess of 60 MeV are temporarily reached. 



panel. Both cases show the characteristic "chirp", up to 
~ 3 ms for the non-spinning and up to ~ 6 ms for the 
tidally locked case, and the subsequent "ringdown" phase 
of the non-axisymmetric central object. The amplitudes of 
the ns^ collisions are displayed in panel two of Fig. [g] Each 
close encounter produces a pronounced GW spike, for ex- 
ample in the {3—1 case, run A, the encounters produce 
spikes at t — 2.3, 8.4 and 12.4 ms. The cases involving black 
holes show substantially longer activity after the first en- 
counter. In the merger cases, run I and J, the episodic mass 
transfer is visible for tens of orbits until the neutron star 
is finally disrupted. The nsbh collision cases are essentially 
"GW-quiescent" (only a small contribution from the close- 
to-breakup rotation of the neutron star core) when the neu- 
tron star is receding from the bh, but produce another GW 
burst at the next close encounter. For the case with the 5 
M0 bh (run E) the longest encountered quiescent phase lasts 



as long as ~ 60 ms. 

The GW amplitude spikes produced by the close encounters 
coincide with peaks in the neutrino luminosities, see below. 
Although the GWs from collisions are comparable in am- 
plitude to those from mergers, the large diversity and the 
lack of a "standard waveform" will make the detection of 
collision signals by current and future ground-based gravi- 
tational wave detectors extremely challenging. 

3.2 Neutrino emission 

We discuss the neutrino properties of one example of each 
encounter class (ns-ns: mergers and collisions; nsbh: mergers 
and collisions) in some detail, an overview over the proper- 
ties of all cases is provided in Tab. |3] Again, we use the 
non-spinning neutron star merger case, run H, as a refer- 
ence point to gauge the other results. Here, the luminosities 
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Figure 3. 3D rendering of the temperature distribution during the grazing impact of two neutron stars (1.3 and 1.4 M©, /3 = 1; run 
A). It is only in the third close encounter (panel 5) that finally a single object forms. In each close encounter a slew of Kelvin-Helmholtz 
vortices forms at the interface between the stars. For display reasons only matter below the orbital plane is shown and the colour bar 
has been restricted to values below 20 MeV. 
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Figure 4. Density cut through the orbital plane at the end of simulation run B (neutron stars with 1.3 and 1.4 M©, /3 = 2). The rotating 
and pulsating central object undergoes a sequence of mass shedding episodes, thereby producing various interacting shocks around the 
central remnant. 



Table 2. Mass distribution: tana is the time at which we analyze the simulation data, mdisk is the resulting 
disk mass, mfb is the mass in fallback material, mesc is the dynamically ejected mass and £^kin the corre- 
sponding kinetic energy. Consistent with our approach in Sect. [33] we capped the numerical velocities in the 
calculation of £^kin,esc at 0.75 c. 



Run tana [ms] mdisk [Mq] mfb [M©] mesc [M©] ^JkinflO^^erg] {vesc) [c] comment 



Collisions 



A 


21.2 


0.27 


0.10 


0.060 


1.15 


0.13 


ns-ns 


B 


9.0 


0.40 


0.10 


0.009 


0.97 


0.22 


ns-ns 


C 


13.2 


0.32 


0.03 


0.030 


3.61 


0.28 


ns-ns 


D 


127.5 


0.24 


0.11 


0.142 


5.70 


0.19 


ns-bh 


E 


143.6 


0.14 


0.04 


0.172 


10.68 


0.24 


ns-bh 


F 


540.3 


0.05 


0.04 


0.134 


8.73 


0.24 


ns-bh 


Mergers 
















G 


20.3 


0.30 


0.06 


0.050 


1.15 


0.15 


ns-ns, corot. 


H 


19.1 


0.25 


0.04 


0.014 


0.23 


0.12 


ns-ns, no spins 


I 


138.7 


0.16/0.22 


0.04 


0.024 


0.61 


0.15 


"disk in disk", inner/both 


J 


139.3 


0.21 


0.03 


0.049 


1.82 


0.18 


ns-bh, no spins 
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Figure 5. Density cut through orbital plane of run E. Panel one (numbering of the panels is from left to right, from up to down) shows 
a snapshot just after the first, panel two after the second and panel five just after the third pericenter passage. Each pericenter passage 
produces a tidal tail. Note that the neutron star core survives even the third pericenter passage. At the end of the simulation it still has 
a mass of ~ 0.1 Mq and moves on a close-to-parabolic orbit away from the black hole. Note that the scales are changing between the 
different snapshots. 
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Figure 6. Gravitational wave amplitudes /i+ (times distance to source; code units = 1.5 km). The upper two panels show neutron 
star neutron star encounters (mergers left, collisions right), the lower two panels show neutron star black hole encounters (mergers left, 
collisions right). 



increase smoothly and peak about 6 ms after contact (at 
t ~ 7.7 ms) with a total of 1.3 x 10^"^ erg/s, see upper panel 
of Fig. |7l The tidally locked case, run G, produces similar 
results, but due to the larger disk mass, see Tab.|2j slightly 
higher luminosities. Since the debris is extremely neutron- 
rich, the neutrino luminosities are dominated by electron 
anti-neutrinos, followed in importance by electron neutri- 
nos and heavy lepton neutrinos (collectively referred to as 
"z^x")? consistent with earlier findings ( Ruffert et al.|p^97 



Rosswog k, Liebendorfer||2003| ). In a recent 2D neutrino- 
hydrodynamics calculation starting from the matter distri- 



bution that resulted from a 3D simulation ( Price & Rosswog 
2006) Dessart et al. ('2009) found that our leakage scheme 
underestimates the heavy lepton neutrino emission, since it 
does not account for the nucleon-nucleon bremsstrahlung 
process. Therefore our heavy lepton neutrino emission re- 
sults are robust lower limits on the true luminosities. In all 
cases apart from run F we find the following hierarchy in 



the mean neutrino energies: {Eux) > (Eu^) > (Eue)- The 
heavy lepton neutrinos are predominantly produced in hot, 
very dense regions and in the exceptional case, run F, the 
densities have already dropped substantially below nuclear 
matter density (pmax < lO"*^*^ gcm""^) at the end of the simu- 
lation when the mean neutrino energies (as given in Tab.[3f 
are measured. 

Even the most gentle ns^ collision with /3 = 1, run A, 
produces a neutrino luminosity that is approximately three 
times larger than the standard double neutron star merger 
case (run H). The merger neutrino lightcurves are rather 
smooth, the collision cases, in contrast, show a much larger 
variability with luminosity changes of up to a factor of two 
on the dynamical time scale of the central object, ~ 1 ms, 
see Fig. [7l middle panel. The more central collisions, run B 
and C, produce neutrino luminosities of about an order of 
magnitude larger than the reference case, see Tab.|3] 
As discussed earlier, we find long-lived, episodic mass trans- 



Mergers vs collisions: the multi-messenger picture 11 



Table 3. Neutrino emission. Iji^ refers to the luminosity of all flavors, the average energies, {E), are measured 
in MeV at time tana, where they are to a very good approximation constant. 



Run tana \ (tana) {Eu^) {E^^) {Eu^) Comment 



Collisions 



A 


16.5 


6.0 X 10^3 


4.3 X 10^3 


10.5 


15.5 


23.5 


ns-ns 


B 


9.0 


3.0 X 10^4 


1.3 X 10^4 


10.1 


16.6 


21.3 


ns-ns 


C 


13.2 


4.1 X 10^^^ 


1.4 X 10^^^ 


11.0 


16.1 


20.8 


ns-ns 


D 


127.5 


3.1 X 10^3 


8.1 X 10^2 


7.1 


10.1 


21.4 


ns-bh 


E 


122.9 


5.3 X 10^3 


5.7 X 10^2 


6.9 


8.5 


14.6 


ns-bh 


F 


540.4 


1.0 X 10^3 


5.6 X 10^1 


9.0 


10.6 


8.9 


ns-bh 


Mergers 
















G 


17.5 


1.8 X 10^3 


1.3 X 10^3 


9.3 


15.5 


27.2 


ns-ns, corot. 


H 


16.0 


1.3 X 10^3 


9.9 X 10^2 


8.0 


14.4 


26.3 


ns-ns, no spins 


I 


138.7 


8.9 X 10^2 


5.7 X 10^2 


6.5 


10.0 


13.8 


ns-bh, no spins 


J 


138.7 


1.0 X 10^3 


6.5 X 10^2 


6.3 


11.0 


14.9 


ns-bh, no spins 



fer in the investigated nsbh cases. During this phase the 
neutrino emission is moderate, but once the neutron star is 
finally disrupted after dozens of orbital revolutions it reaches 
peak values of up to 10^^ erg/s. All investigated nsbh colli- 
sions robustly produce hot and massive accretion disks with 
neutrino peak luminosities of at least 75% of the standard 
merger case. We display the neutrino luminosities for the 
case shown in Fig. [5] run E, in the last panel of Fig. [7| All 
pericenter passages are visible as a peak in the neutrino lu- 
minosities. Passages two and three thereby substantially in- 
crease the disk mass, the subsequent disk consumption each 
time produces an enhanced neutrino luminosity. In the case 
of run E, the enhanced emission due to a pericenter passage 
is therefore double peaked with a first, short peak resulting 
from the neutron star remnant impact onto the disk (be- 
tween panels four and five in Fig. [5]), and a second, broader 
peak driven by the viscous consumption of the re-filled ac- 
cretion disk. In summary, each pericenter passage produces 
both a gravitational wave and at least one neutrino emission 
peak. 



3.3 Fallback and external feeding of accretion 
disks 

Swift observations have revealed that both short and long 
bursts can show late time X-ray flaring activity (Burrows' 
& et ah, 2005; Nousek & et al. 2006). Some SWIFT sGRBs 



are accompanied by a phase of extended X-ray emission that 
lasts between - 10 and 100 s (Norris Gehrels,2QQ8, : ,Perley, 
et al.|[2009t and whose fluence can exceed that of the GRB 
itself. Such time scales are substantially longer than the dy- 
namical/viscous time scales, i.e. ~ 1 ms/0.05 s, that are ex- 
pected for a merger remnant. In this context, it was pointed 
out that compact binary mergers possess in addition a much 
longer time scale due to nearly unbound material that will 
finally fall back to the central remnant ( Rosswog|2007a Lee 



Ramirez-Ruiz||2007| |Faber et al^|2Q06 |). The mass in the 



fallback of a merger however is substantially lower than that 
of the accretion that is supposed to launch the GRB. There- 
fore it is at least not obvious how fallback would produce 
extended emission with a larger fluence than the sGRB. 



Clearly, the accretion of the fallback material is a rather 
complicated phenomenon and much of what will flnally be- 
come visible as electromagnetic radiation will depend on the 
evolution of the inner remnant disk as it interacts with the 
fallback material being delivered at super-Eddington rates. 
It has been argued ( Rossi Begelmaii||2009 ) that -rather 
than being flnally swallowed by the likely emerging central 
black hole- the fallback material could form an extended hot 
envelope around the remnant once neutrino cooling ceases 
to be an eflicient cooling agent, but the matter is still opaque 
to photons. After about one week a soft X-ray signal near 
the Eddington luminosity should become visible. 
In addition, while the neutron star debris is initially receding 
from the central remnant it will be subject to heating from 
radioactivity ( Li Paczyhski|1998| Rosswog|2005 Metzger 
et al.|[20T0 Roberts et al.||20lT ) and it has been proposed 
that this radioactive energy release might regulate the fl- 



nal fallback rate (Metzger et al. 2010| ). At the same time, 
viscous effects drain the disk via accretion onto the central 
object and the disk expands due to outward angular momen- 
tum transport. Winds are expected as a result of viscous 



heating and/or neutrino heating (Lee et al. 2005 Dessart 



et al.|2009J which further complicates the situation. At late 



times, recombinations of nucleons into a-particles are ex- 
pected to unbind a substantial fraction of the late-time disk 



( |Beloborodov|2008| [Lee et al.||2009| [Metzger et al.|2009 ). 

Clearly, this interesting, but complicated evolution deserves 
further detailed studies, but reliable predictions of the elec- 
tromagnetic display require a reliable knowledge of the prop- 
erties in the emerging photosphere. The treatment of these 
effects is beyond the scope of the current work and we re- 
strict ourselves here to a simple semi-analytical model to 
estimate the time scales on which fallback material is de- 
livered to central object while ignoring for now additional 
complicating effects such emerging winds. 
We assume that for most of its trajectory, the fallback ma- 
terial is reasonably well described as being ballistic. Un- 
der this assumption we can calculate quantities like eccen- 
tricity, circularization radii etc. as described in detail in 
Rosswog' ('2007a) and we can in particular calculate the 
time, tfb(i?circ), it takes each matter portion to return to 
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Figure 8. The dependence of the matter energy accretion rate on 
the Shakura-Sunyaev viscosity parameter oigs is illustrated for the 
case of a merger between a ns and a 5 bh (run I). Increasing 
the viscosity time scale (= reducing the parameter ass) smoothly 
spreads out the accretion to longer time scales. 



its circularization radius, i^circ- Under the simplifying as- 
sumption that matter quickly settles into a disk, we assume 
that from this point of time, tfb,a, each fluid parcel (SPH 
particle) a is accreted at a constant rate jria = ma/tvisc,a 
on the viscous time scale of a thick accretion disk (with 

H ^ R), tvisca = (ass^K(i?circ,a))~\ whcrC Qss IS the 



Shakura-Sunyaev viscosity parameter (Shakura Sunyaev 
1973). The total mass accretion rate is then given as a sum 



over particle contributions 
M(t) - ^ma(t). 



(10) 
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Figure 7. Examples of the emerging neutrino luminosities: the 
"standard" merger case (run H, upper panel), the ns"^ collision 
with /3 = 1 (run A, middle), and the nsbh collision with /3 = 1 
and mbh = 5 (run E, bottom). 



where 



^a(t) 



J 'vise, a 

I 0. 



tfb,a 

else 



< t ^ tfb,a + Tv: 



(11) 



This simple model extends to original one ( |Rosswog|20Q7a| ) 
by accounting also for the viscous dissipation time scale. Ad- 
mittedly, this model is very simple and the complicated but 
important topic fallback accretion in the aftermath of a com- 
pact binary merger deserves more efforts in future work. 
We illustrate the effect of the viscosity parameter a^s in 
Fig. [8] for the case of a merger between a ns and a 5 
bh (run I). We show the resulting curves Mc^it) for cases 
where either viscous time scales are ignored (i.e. a^s = oo) 
or varied in a plausible range (o^ss = 0.1, 0.01 and 0.001). 

In Fig. [9] we show the values of Mc^ from fallback ma- 
terial as predicted from our model. For all shown curves 
a value of ass = 0.01 was adopted which is realistic, but 
maybe slightly on the low side. In our standard merger 
case, run H, the fallback luminosity peaks around 0.1 s with 
^ 3 X 10^^ erg/s. The somewhat academic corotating case, 
run G, shows a similar peak fallback rate which is nearly 
constant for the first ^ 0.15 s. The /3 = 1 ns^ collision case, 
run A, produces a roughly constant fallback rate until 0.15 
s, before it enters the 5/3-powerlaw phase that is expected 
for dM/de ^ const ( |Rees||1988| |Phiniiey||1989i ). During this 
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Figure 9. Mc^ from fallback material as predicted from the simple model (ballistic fallback + disk accretion) described in the main 
text. Left: neutron star - neutron star cases, right: neutron star - black hole cases. In each panel both merger and collision cases are 
shown. For comparison we also show the power law for "standard" fallback (Rees 1988 Phinney 1989) that assumes a constant mass per 
energy, dm/de = const, and results in a time dependence oc t~^^^. For the figures a viscosity value of ass = 0-01 has been adopted, see 
text for details. 



phase the fallback luminosity is about a factor of four larger 
than in the standard neutron star merger case. The other 
ns^ collision cases peak earlier and their fallback rates are 
lower than for canonical ns^ mergers. The early phases of 
the fallback curves should generally be interpreted with a 
grain of salt since they are sensitive to the exact location of 
the disk radius and the assumption of ballistic motion is not 
necessarily well justified. At later times (t > 0.2 s) all cases 
are fairly well described by the M oc t~^^^ powerlaw which 
is shown in the figure as a dashed line. 

The nsbh cases are shown in the right panel. The mergers, 
run I and J, are qualitatively similar to the ns^ merger cases, 
after a "plateau phase" up to ^ 0.3 s, the fallback luminos- 
ity drops off compatible with M oc t~^^^ (the lack of late 
time data points for nsbh mergers is due to the substantially 
lower numerical resolution for these cases). The nsbh colli- 
sions show less clear initial plateau phases and only settle 
into the 5/3-powerlaw phase after a few seconds. 
The collision cases yield typically about a factor of two more 
fallback material than the mergers, see Tab. [2] This mass 
is located in tidal tails (one per pericenter passage) which 
can be considered as external mass reservoirs that keep re- 
plenishing the central accretion disks at their fallback rates. 
With an available energy of ^ 2 x lO^^erg (e/0.1) (mfb/0.1 
M0) the fallback material is in principle able to produce en- 
ergetic post-merger activity and, as evident from Fig. [5] the 
fallback time scales are orders of magnitude longer than the 
dynamical time scales. But as discussed above, the real situ- 
ation is possibly too complicated to be adequately captured 
by this simple fallback + disk model. Therefore it remains 
an open question whether fallback is the trigger of the late 
activity that is observed in some sGRBs or not. 



3.4 Ejecta 

All investigated cases dynamically eject a substantial 
amount of neutron-rich material, see Tab. These dy- 
namical ejecta are complemented further by neutrino-driven 
winds (Dessart et al.'2009 ) and late-time evaporations of the 
accretion disks ( Beloborodov|2008j Lee et al. 2009 i Metzger] 
et al.||2009 ). The amount of ejecta that we find in the pre- 
sented simulations is very robust with respect to changing 
the numerical resolution, it is essentially unchanged even if 
we reduce the particle number by an order of magnitude. 
The amount of dynamical ejecta may change, however, with 
a different treatment of gravitational wave backreaction (re- 
duction) , general relativity (probably reduction) or a differ- 
ent nuclear equation of state (depending on EOS stiffness 



either reduction or enhancement). In earlier work ( Ross wog| 
et al.||2000 ) we had seen the tendency of stiff nuclear equa- 
tions of state to eject more material than softer ones. The 
recent discovery of a neutron star close to 2.0 Mq (Demor- 



est et al.||2010|) suggests that the true neutron star EOS is 



indeed very stiff, therefore we consider the Shen et al. EOS 
as a good choice. 

The amount of dynamically ejected material is summarized 
in Tab. [2] The most likely case, run H, ejects 1.4% of a solar 
mass with a kinetic energy of 2 x 10^^ erg, the corotating 
neutron star merger case ejects three times as much mass. 
The collision cases eject substantially more, about four times 
as much for the most likely ns^ collision, run A, and up to 10 
times as much for nsbh collision, run E. The kinetic energies 
in the ejecta can reach up to 10^^ erg (run E). 
The ejected material is extremely neutron-rich and origi- 
nally resides in the crust and the outer core of the neu- 



^ An overview over the ejecta masses for more than 30 different 
cases can be found in Rosswog ( 2012| ). 
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Figure 10. 3D rendering of the Ye distribution for the generic neutron star merger case (1.3 M©, 1.4 M©, no spins; run H). Only matter 
below the orbital plane is shown. The bulk of the dynamically ejected matter has Ye-values as low as ~ 0.03. 



tron star. We display a 3D rendering of the Ye distribution 
for our standard neutron star merger case in Fig. (val- 
ues are shown at a given artificial optical depth, see [Price] 



(2007)). The fractional Ye distribution within the eject a of 
the different cases is displayed in Fig. [TT] Both neutron star 
merger cases, runs G and H, upper left panel, eject material 
with a pronounced peak near Ye ~ 0.03 with (marginally re- 
solved) higher Ye-contributions from the neutron star crust, 
see Fig. [ij For the most likely ns^ collision, run A, the Ye 
distribution looks similar (upper left panel). The more ex- 
treme -and less likely- more central collision cases {/3 — 2 
and 5) produce strong shocks with very large temperatures 
{Tpeak ~ 80 MeV) . In these cases positron captures substan- 
tially increase the electron fraction. In both these cases the 
Ye distribution has a broad peak near Ye ^ 0.2 (see upper 
right panel). 

The nsbh cases show a similar tendency: where the eject a are 
predominantly due to tidal torques they still possess their 



original, very low Ye. If, in contrast, they suffered strong 
heating in shocks, the Ye distribution is shifted to larger Ye 
values. The nsbh cases experience a different dynamical evo- 
lution. In the long episodic mass transfer phase, see Fig. [6] 
they are continuously tidally heated. In the later stages of 
mass transfer a disk starts to emerge and there is a contin- 
uous hydrodynamic interaction between this growing disk 
and the ns remnant core. Due to this heating history the Ye 
distribution is different: it still shows a peak near Ye — 0.03 
but also extends smoothly to values of ^ 0.4. The collision 
cases also show somewhat increased Ye values (^ 0.1) that 
are due to the various heating phases in the different peri- 
center passages. 

The distribution of asymptotic velocities within the dynamic 
ejecta is shown in Fig. 12 The distributions in ns^ merger 



cases (upper left panel) show peak values around 0.15 c and 
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Figure 11. Ye distribution in the dynamical eject a binned by mass. The upper row shows double neutron star encounters (mergers left, 
collisions right), the lower row shows neutron star black hole encounters (mergers left, collisions right). 



hardly any material faster than 0.4 The nsbh merger 
cases show very similar velocity distributions. The collision 
cases, in contrast, show a much broader distribution of ve- 
locities with small mass fractions reaching beyond 0.7 c, see 
the right column of Fig. |12[ The bulk of ejecta in the graz- 
ing ns^ collision has velocities similar to the merger cases 
[r^ 0.15 c), but smaller amounts of matter can reach sub- 
stantially higher velocities. The more central collisions show 
high velocity tails reaching to close to the speed of light. 
Such close-to-c velocities in the central collision cases (/3 = 2 
and 5) are likely an artifact of our essentially Newtonian 
simulations. The nsbh collisions show a very homogeneous 
behavior with the bulk of matter reaching about 0.3 c and 



^ We chose the lower limit of the ordinate depending on the nu- 
merical resolution. Since the nsbh merger cases have a lower res- 
olution this value differs from the nsns cases. 



peak velocities out to 0.8 c. The impact of these properties 
on the electromagnetic signal is further discussed in Sec.|3.5| 



3.5 Electromagnetic emission 

S. 5. 1 Radio remnants from ejecta-ISM interactions 

At large radii the sub and mildly relativistic outflow from 
mergers/collisions is decelerated by the external medium, 
driving a fast shock into it. Shocks with similar velocities 
are seen in the late phases of GRBs and in early phases of 
some supernovae, and are known to produce a bright radio 
emission. This radio emission is explained by the large frac- 
tion of the internal energy behind the shock, ^ 10%, that is 
channeled into accelerated electrons and a similar amount 
of energy that goes into magnetic fields. In |Nakar &: Piran] 
( |2011| ) and in Paper II (| Piran et aL]|2012 ) we describe the 
calculation of the radio emission that follows compact bi- 
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Figure 12. Distribution of the asymptotic velocity within the dynamical ejecta. Upper row ns"^ encounters, mergers left, collisions right. 
Lower row, nsbh encounters, mergers left, collisions right. 



nary mergers. Here we use the same method to compare the 
radio remnants of mergers and collisions. 

In Figs. [13] and [M] we present the radio signal that re- 
sults from the dynamical ejecta being slowed down by the 
external medium. The calculation is based on the velocity 
profiles presented in Fig. [12] and it approximates the in- 
teraction as a spherical blast wave in an ambient medium 
with a constant density, n. Behind the shock constant frac- 
tions of the internal energy, te = 0.1 and cb = 0.1, are de- 
posited in relativistic electrons and in magnetic field, where 
the electrons are accelerated to a power-law with an index 
p = 2.5. The light curves are calculated following the pro- 
cedure described in Paper II ( |Piran et al.||2012 ). Since the 
simulations are Newtonian we do not trust velocities that are 
close to the speed of light and we therefore conservatively 
restrict the simulation velocities to a maximum of 0.75 c 
when using them to calculate the radio emission. Note that 
the simulations presented here do not include all possible 



sources of outflows from such events, and especially sources 
of mildly relativistic and relativistic ejecta (e.g., various out- 
flow sources from near the event horizon of the accreting 
black hole, such as Blandford-Znajek fBlandford & Znajek 



1977)) are not accounted for. Thus, the true remnants may 



be brighter than our prediction, especially on short time 
scales (weeks-months) where the dissipation of the energy 
in the fastest moving ejecta is most efficient. 
The radio lightcurves of ns^ and nsbh mergers/collisions, 
see Figs. [13] and [M] are calculated at two frequencies, 1.4 
GHz and 150 MHz, and for two values of external densities 
n — 1 cm~^ and n = 0.1 cm~^. The flux normalization is 
for events at a distance of 10^^ cm, roughly at the detec- 
tion horizon for ns^ mergers by advanced LIGO and Virgo. 
The main difference between collisions and mergers is that 
collisions produce brighter radio remnants that rise faster. 
Among the collisions the peak flux increases with the pene- 
tration factor 13. In general the flux from collision remnants 
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Figure 13. Radio lightcurves of ns^ mergers and collisions at 150 MHz and 1.4 GHz for two different densities of the external medium 
(1 cm""^, left and 0.1 cm~^, right). The microphysics parameters are ee = = 0.1 and p = 2.5, the flux normalization is for events at 
a distance of 10^^ cm, roughly the detection horizon for advanced LIGO and VIRGO. 
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Figure 14. Same as Fig. [13] but for the encounters involving a black hole. 



rises fast, within less then a yeai[^ and remains bright (or 
even continues to rise, but more slowly) for several years 
for ns^ collisions and for decades in the nsbh collision case. 
Merger remnants, in contrast, are fainter and they rise over 



^ The rise time is set in this case by the artiflcial cut-off (0.75 c) 
that we took in the outflow velocity. A signiflcant outflow with 
higher velocity results in a shorter rise time. 



time scales ranging from a year to a decade and remain 
bright for a decade or more. The reason for these differ- 
ences is that collisions are producing more energetic and 
faster outflows. The larger energy and the faster ejecta in- 
crease the remnants' brightness while the faster ejecta are 
the dominant factor in the shorter rise time, see the detailed 
discussion in Paper II (Piran et al.|201 2) of the effects of the 
outflow's energy and velocity on the radio signal. Interest- 
ingly, a lower external density reduces the flux and slows its 
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evolution, as long as the observed frequency is above the syn- 
chrotron self-absorption at all time. This is the case at 1.4 
GHz and therefore, observing at this frequency alone, one 
cannot distinguish between a ns^ merger at n = 1 cm""^ and 
a collision (with ^ — 1) that take s pl ace within a medium 



of density n — 0.1 cm (see Fig. 13). The key for remov- 



ing the degeneracy is to observe the remnant also at 150 
MHz, where the effect of the external density on the self- 
absorption frequency leaves its mark on the lightcurve and 
makes the two scenarios (merger in dense environment vs. 
collision in sparse environment) distinguishable. 

The major uncertainty concerning the detect ability of 
collisions (the detectability of mergers is discussed in Pa- 
per H ( ,Piran et al.||2012" )) is due to the external density in 
globular cluster cores. While globular clusters usually reside 
in very low density environments, the gas in their cores is 
expected to be dominated by the mass loss from post main 
sequence cluster members. The only globular cluster where 
gas is clearly detected is 47 Tuc, where the d ensity of ion- 
ized gas is of order 0.1 cm~^ (Freire et al. 2001 ). This density 
corresponds to a total mass of 0.1 Mq in the cluster core. 



van Loon et al. ( 2006 J reports on a tentative detection of 
0.3 Mq of neutral hydrogen in M15 and possible detections 
of similar amount of neutral hydrogen in two other clusters 
(and upper limits of ^ Mq in two other). Therefore, it is 
realistic to expect a circum collision density of 0.1 cm""^, or 
even larger, nsbh cases are expected to dominate the radio 
remnant population from collisions, since their volumetric 
rate is higher ( Lee et al.|2010 ) and they are about ten times 
brighter. Since these remnants remain bright for decades the 
best opportunity to identify a remnant is to catch it dur- 
ing its fast rise phase, which takes about a year. Assuming 
the upper limit for nsbh collisions (6 x 10~^ per Milky way 
galaxy per year, see below) and taking a detection distance 
of 1 Gpc (as expected for a sub-mJy survey at 1.4 GHz) the 
number of remnants which are younger than 1 year, over the 
whole sky, is < 200. 



3.5.2 Transients from radioactive decays, "macronovae^^ 

As discussed in the previous section, the ejected material is 
extremely neutron-rich and rapidly expanding. Under such 



conditions rapid neutron capture is hard to avoid (Hoffman 
et al.| |1996'; "Freiburghaus et al.|[T999l [Roberts et al.||2011| ). 

The r-process itself occurs on a very short time scale, but the 
freshly synthesized elements subsequently undergo nuclear 
fission, alpha- and beta-decay which occur on much longer 
time scales. 

The supernova- like emission powered by this radioctive 
decay of this expanding material was first suggested by_0 



& Paczyhski (1998) and discussed later by Kulkarni (2005) 
and Metzger et al. (2010). Here we use the ejecta velocity 



profile provided by our simulations to calculate the observed 



light curve following the formalism of Piran et al. ( 2012| - 
paper H. 

The resulting macronova lightcurves are shown in Fig.|15|for 
both ns^ (left) and nsbh encounters (right). The canonical 
merger case (left, black) peaks after 0.4 days with ^ 5 x 10^^ 
erg/s. The collision cases show a spread in peak times from 
0.1 to 1 days and peak luminosities up to 10^^ erg/s. The 
macronova resulting from the merger of the 1.4 M© ns with a 
5 M© bh (right, black) is similar to the canonical ns^ merger 



case, the 10 M© nsbh merger peaks at about 0.8 days with 
10^^ erg/s. The nsbh collision cases all peak beyond 1 day 
with - 2 X 10^^ erg/s. 



4 SUMMARY AND DISCUSSION 

Dynamical collisions between two neutron stars and a 
neutron star and a stellar mass black hole are naturally 
expected in stellar systems with large number densities 
such as globular clusters. A recent study ( Lee et aL]|2010 ) 
found that dynamical collisions between two neutron stars 
and a neutron star and a stellar mass black hole occurs 
more frequently than previously estimated and that such 
collisions could contribute to the observed sGRB rate. 
In a large set of simulations we have explored the multi- 
messenger signatures that are produced by both dynamical 
collisions and the more conventional gravitational wave 
driven compact binary mergers. These simulations use 
Newtonian gravity, but benefit from the use of a nuclear 
equation of state, a multi-flavour neutrino treatment and a 
numerical resolution (up to 10^ SPH particles) that exceeds 
existing studies by far. Nevertheless, this study will need to 
be updated once relativistic simulations with microphysics 
become available. 



Dynamics 

Generally, collisions show a larger variety in all of their 
properties than binary mergers, mainly due to the lack 
of strong constraints on their mass ratio and the impact 
parameter. Typical velocity dispersions in globular clusters 
are orders of magnitude smaller than the velocities the 
compact objects reach due to their mutual gravitational 
attraction, therefore initially parabolic orbits are good 
approximations. A compact binary system is already 
strongly bound at the onset of the dynamical merger phase, 
collisions, in contrast, possess a total orbital energy close 
to zero. To finally form a single remnant they have to get 
rid of energy and angular momentum by gravitational wave 
emission and mass shedding episodes. Colliding systems 
therefore undergo several close encounters before they can 
finally merge into a single central object surrounded by 
debris. Each close encounter launches a tidal tail and the 
neutron star(s) is/are spun up to rotation frequencies close 
to breakup. The final remnant is then, like in the merger 
case, a supermassive neutron star/black hole surrounded by 
a disk, but the disk is externally fed by tidal tails (one per 
close encounter). A good fraction of the mass is still bound 
to the central remnant and will fall back on times scales 
substantially exceeding those of the central engine. 

GRB engines 

We find that collisions are at least as promising as GRB 
engines as gravitational wave driven binary mergers. For 
example, ns^ collisions naturally produce Kelvin-Helmholtz- 
unstable shear layers in each close encounter, see Fig. [TgJ 
which have been found in previous to amplify initial neutron 
star magnetic fields ( Price Rosswog|2006 Anderson et al. 



2008 Obergaulinger et al. 2010) and they robustly form 
massive accretion disks between ~ 0.05 to ~ 0.4 M©. 
Their neutrino luminosities are at least comparable to the 
standard neutron star merger case, ^ 10^^ erg/s, in central 
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Figure 15. Radioactive decays within the dynamical ejecta produce transient events ( "macronovae" ) in the optical-UV spectral range 
that share several properties with supernovae, but evolve on substantially shorter time scales. The lightcurves for ns"^ encounters (mergers 
and collisions) are shown in the left panel, those of nsbh encounters in the right. See main text for details. 




Figure 16. Kelvin-Helmholtz unstable shear interfaces are a natural consequence of ns^ collision (illustrated here in 3D rendering of the 
temperature distribution; only matter below orbital plane is shown). Left: grazing inpact (13 = 1, run A), right: more central collision 
(^ = 2, run B), colour bars restricted to T < 20/24 MeV to enhance visibility. 



ns^ collisions this value can be exceeded by more than an 
order of magnitude, see Tab. [3] The engine dynamics leaves 
a clear imprint in the neutrino luminosity. The oscillations 
of the freshly formed central remnant of a ns^ collision 
show up as ms oscillations in the neutrino luminosities, 
see Fig. |7| middle panel. In nsbh collisions the core of the 
neutron star can survive several close encounters. Each of 
them "re-fills" the accretion disk and thus produces apart 
from a GW burst also peaks in the neutrino luminosities. 
We find examples of double peaks where the first occurs 



when the ns core impacts on the disk and the second results 
from the consumption of the "refilled" disk, see for example 
run E (ns of 1.3 M© and bh of 5 M©), Fig.[7| last panel. 
If neutrino annihilation should be the main driver for 
launching relativist ic outflows, collisions would be at least 
as promising as binary mergers. One caveat, though, that 
applies to both mergers and collisions, is that neutrinos 
can also deposit part of their energy (either via absorption 
or annihilation) in the remnant matter and thus drive 
strong baryonic winds ( Dessart et al.|2()09 ). The latter work 
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found that for at least as long as the central neutron star 
has not yet collapsed into a black hole a strong baryonic 
wind is blown along the original rotation axis which may 
prevent relativist ic outflow from forming. This threat may 
be specific for cases where a central neutron star survives 
for many dynamical time scales. In such cases neutrinos 
deposit energy in the outer neutron star layers and thus 
blow off baryons, see Fig. 10 in Dessart et al. (2009). To 
which extent this effect also occurs in nsbh encounters 
is currently not known. This interesting topic certainly 
deserves more detailed work in the future. 

Fallback and late-time activity 

A particularly attractive feature of the collision hypothesis 
is the possiblity to launch large amounts of matter into 
eccentric "fallback" orbits. Collisions can deliver more 
than twice the fallback mass of the standard neutron 
star merger case, see Tab. |2] Our fallback plus thick disk 
accretion model shows time scales in excess of 10^ s, but 
it is too simple to account for possible flares. Generally, 
Mc? has dropped below lO'^^ erg/s at times beyond 10^ 
s. An intriguing possibility, though, is the survival of the 
neutron star core in nsbh collisions. In all investigated cases 
it survives at least the first encounter, in one calculation 
even three of them and is finally ejected with 0.1 in 
a close-to-parabolic orbit. The interaction with the debris 
from previous passages may actually brake the core enough 
to fall back to the bh. This would trigger an impulsive 
accretion event with ^ 10^^ erg at possibly very late 
times. Such scenarios may be responsible for the observed 
late- time activity in some cases, but they are probably 
too rare (see below) to account for extended emission in a 
substantial fraction of sGRBs. 

Detectability of the neutrino signal 

Most of the energy of both mergers and collisions escapes as 
^15 MeV neutrinos. The neutrino signature is comparable 
in magnitude and duration to the signals of a typical 
core collapse supernova, though dominated by electron 
antineutrinos. Current facilities can only detect individual 
neutrino events from within the galaxy (or at most events 
within the local group) . Even if more sensitive detectors are 
built in the distant future it will be difficult to distinguish 
the very rare merger or collision events among the much 
more frequent signals from core collapse SNe. So the 
neutrinos from compact object encounters may just make 
an individually undetectable ^ 0.1% contribution to the 
diffuse neutrino background. 

Detectability of the gravitational wave signal 
In both mergers and collisions gravitational radiation 
carries also a significant fraction of the energy, second 
only to the neutrino signal. The quasi-regular mergers' 
chirping signal can be detected up to distances of a few 
hundred Mpc, they are the prime targets for gravitational 
radiation detectors like advanced LIGO and Virgo. Since 
gravitational waves from eccentric binary systems efficiently 
radiate angular momentum relative to energy, mergers of 
primordial binaries are expected to occur at practically zero 
eccentricity. Binaries that have formed dynamically, say in 
nuclear or globular clusters, however, form at small orbital 
separations but with large eccentricities and they may not 



have enough time to circularize up to merger. Thus, their 
close-to-encounter orbital dynamics and gravitational wave 
signal may differ substantially. Initially they produce a 
series of well- separated, repeated GW bursts that continues 
for minutes to days. This sequence of bursts gradually 
transforms into the powerful chirp inspiral signal of an 
eccentric binary system. Kocsis et al. (2006) found the 



signal-to-noise ratio for GWs produced in single parabolic 
passages to be significant only for rather deep encounters 
where the initial pericentre separation is rp,o < 6M, M 
being the total mass (G = c = 1). O'Leary et al. (2009) 



estimate that nsbh mergers that result from tidal capture 
binaries in galactic nuclei could be detectable by advanced 
LIGO type facilities at rates of ~ 1 per year, i.e. they could 
contribute of order 1 % to the overall detection rate. Since 
a significant fraction of tidal capture binaries is expected 
to merge at non-negligible eccentricities their properties 
are bracketed between those of the quasi-circular binary 
mergers and the 13 ^ 1 collisions that were discussed in 
this paper. A first attempt to model such high eccentricity 
mergers has been undertaken by [East et al. (2012), further 
modeling is left for future efforts. 



Rate constraints from ejected matter 

It had long been suggested that the cold decompression of 
neutron star material could produce a substantial contribu- 



tion to the r-rprocess inventory of the Unive rse (Lattimer 



Schramm] [T974| |T976l [Eichler et al.||1989| [Freiburghaus 
et al.||1999||Rosswog et al.||1999| ). Core-collapse super novae 
are often considered as the "standard" r-process production 
site, but recent studies suggest that they are seriously chal- 
lended in producing the whole observed r-process pattern 
( Roberts et al.|[20To Fischer et al.||20l"o{ [Arcones Janka 



2011) unless very specific SN events are invoked (Winteler 



et al. [2012 ). These findings make an alternative/additional 
r-process production channel more welcome than ever. 
Our reference case, a neutron star merger with masses of 
1.4 M0 (run H), ejects 0.014 Mq of extremely neutron- 
rich (y'e'^ 0.03) matter, mainly from the outer core with 
small contaminations from the more proton-rich crust. 
Grazing ns^ collisions (run A) eject material 
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see Fig. 

with very similar properties, while more central collisions 
yield substantially higher temperatures/e^-capture rates 
and therefore larger Ye values {Ye^ 0.2), though still very 
low in comparison to core-collapse SNe. If we assume an 
average Galactic r-process production rate of 10~^ Mq per 



year ( |Qian 2000) we can place robust upper limits on the 



admissible encounter rates under the assumption that the 
ejecta are made entirely of r-process nuclei. Our recent 
study ( Korobkin et al.|2012 ) finds indeed that the dynamic 
ejecta of ns and nsbh mergers are excellent candidates for 
the production sites of the robust heavy r-process patter 
with A larger than ^ 120. If neutron star mergers alone were 
responsible for all r-process material this would constrain 
their galactic rate to 7 x 10~^ M© yr~^, right in the centre 
of the 95% confidence i nterval derived from the observed 
ns^ binary distribution (Kalogera et al.|[2004), see Fig 
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Within the given uncertainties they would be perfectly 
consistent with delivering a substantial contribution to the 
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In analogy with supernova remnants, the matter that is 
ejected at sub and mildly relativistic velocities produces 
a longer lasting radio flare (Nakar & Piran 2011). The 
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Figure 17. Nucleosynthetic constraints on the event rates. The 
double neutron star merger rate that is required in order to pro- 
duce all the r-process material (using the average rate of 10~^ 
M0 yr"-*^, dashed line, Qian 2000) is marked by the filled circle. 
This rate is promisingly consistent with the rate derived from ob- 
servations by Kalogera et al. (2004) shown as green band. If we 
make the (extreme) assumption that only ns^ and nsbh collisions 
produce all the r-process and their relative ratio is 1:5 (Lee et al. 
2010) we can place upper limits on their occurence rates. These 
limits are marked as triangle (ns^) and square (nsbh). The rates 
realized in nature might be substantially below these values. 



galactic r-process matteijj The collision cases with their 
larger ejection yields, however, are seriously constrained 
in their allowed occurrence rates. If we make the extreme 
assumption that ns^ and nsbh collisions produce all the 
r-process and use a nsbh rate, i^coii,nsbh, that is five time s 
larger than the ns^ collision rate i^coii,ns2 ( Lee et a'L||2010 ), 
we find upper limits of i?coii,ns2 < 1.2 x 10~ yr~ and 

,nsbh 

< 6 X 10"^ yr"^ Therefore the total (ns^+nsbh) 
collision rate must be i?coii < 7.2 x 10~^ yr~^, or about 10% 
of the double neutron star merger rate. Given our extreme 
assumption that no other astrophysical event produces 
r-process material the true rate is likely substantially below 
this estimate. 



Macronovae 

We estimate the properties of the electromagnetic tran- 
sients due to radioactive decays for spherically symmetric 
outflows. The standard neutron star merger case peaks after 
~ 0.4 days with a luminosity of 5 x lO'^^ erg/s. The other 
merger cases (ns^ and nsbh) and the ns^ collisions peak 
with higher luminosities, both within a factor of about two 
in comparison to the standard merger case. The neutron 
star black hole collisions form a distinct group: they all 
peak beyond 1 day and are substantially brighter that the 
standard merger case (^-^ 2 x 10^^ erg/s). 

Radio signal from ejecta-ISM interactions 



^ It remains to be investigated, though, to which extent this is 
consistent with galactic chemical evolution, see e.g. ,Argast et al.| 
120041) 



flare is the most robust electromagnetic counterpart that 
is expected from compact mergers or collisions. It depends 
only on the total energy ejected in mildly or subrelativitic 
energies and on the density of the external matter. A flare 
from a canonical ns^ merger would peak on a time scale 
of a year with a peak observed flux of ^ 0.1 mJy at 1.4 
GHz from a source from a distance of lO^^cm, roughly 
the detection horizon of the gravitational radiation signal 
by the advanced LIGO/Virgo. The signals are longer 
and weaker if the external density is lower, ns^ collisions 
or bhns collisions eject much more material and with a 
higher velocity and hence the corresponding flares would by 
stronger and longer. In addition the rise time of these flares 
will be faster (because of the higher velocities) and hence 
it will be easier to detect them. In general the rising phase 
of the flares depends critically on the amount of fastest 
matter ejected. In these Newtonian calculations we had 
capped the maximal velocity at 0.75 c hence our results on 
the early emission underestimate the true emission. This 
part of the relativistic ejecta that we miss in the simulation 
may produce a brighter and a faster rising signal even if its 
energy is ten times smaller. 
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